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Abstract 

We present a new method for imposing a realistic equation of state in anisotropic hydrodynamics. 
The method relies on the introduction of a single finite-temperature quasiparticle mass which is 
fit to lattice data. By taking moments of the Boltzmann equation, we obtain a set of coupled 
partial differential equations which can be used to describe the 3-|-ld spacetime evolution of an 
anisotropic relativistic system. We then specialize to the case of a 0+ld system undergoing boost- 
invariant Bjorken expansion and subject to the relaxation-time approximation collisional kernel. 
Using this setup, we compare results obtained using the new quasiparticle equation of state method 
with those obtained using the standard method for imposing the equation of state in anisotropic 
hydrodynamics. We demonstrate that the temperature evolution obtained using the two methods 
is nearly identical and that there are only small differences in the pressure anisotropy. However, 
we find that there are significant differences in the evolution of the bulk pressure correction. 
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I. INTRODUCTION 


Ultrarelativistic heavy-ion collision experiments using the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory and the Large Hadron Collider (LHC) at CERN 
allow researchers to study the behavior of matter subject to extreme conditions. In these 
experiments, high-energy collisions of nuclei are used to heat a tiny volume of matter up 
to temperatures that exceed the critical temperature (Tc ~ 160 MeV) necessary to create a 
super-hot deconhned and chirally-symmetric phase, called the quark-gluon plasma (QGP). 
The study of this strongly interacting phase near and above the critical temperature is of 
fundamental interest. One can gain some insight into the physics of the QGP using per¬ 
turbation theory since the asymptotic freedom of quantum chromodynamics (QCD) ensures 
that, for the high temperatures, T 3> Aqcd, the QGP can be thought of as a weakly- 
coupled many-body system. In this regime, perturbative methods, such as hard thermal 
loop (HTL) resummation, can be used [1-6].^ In the HTL framework, the quarks and glu¬ 
ons can be thought of as quasiparticles having temperature-dependent (thermal) masses 
with rrig^g^g ~ gT, where g is the strong coupling. 

Such a picture provides motivation to try to model the QGP as a gas of massive quasipar¬ 
ticles for the purposes of obtaining self-consistent hydrodynamic equations. However, per¬ 
turbation theory needs to be supplemented since, for temperatures T < 2Tc, hrst-principles 
perturbative calculations based on deconhned quarks and gluons break down. In order 
to proceed, one can use non-perturbative lattice calculations of QCD thermodynamics to 
determine information about the necessary quasiparticle mass(es). In practice, one can 
perform this procedure at all temperatures and determine a non-perturbative temperature- 
dependent quasiparticle mass, m(T). Once m{T) is determined, one can use this to enforce 
the target equation of state (EoS) in an effective kinetic held theory framework. One compli¬ 
cation is that, in order to guarantee thermodynamic consistency in equilibrium and related 
out-of-equilibrium constraints, it is necessary to introduce a background (vacuum energy) 
contribution to the energy-momentum tensor [9-11]. The resulting EoS, together with a self- 
consistent non-equilibrium energy-momentum tensor and modihed Boltzmann equation, can 
be used to derive relativistic hydrodynamic equations for such a quasiparticle gas. 

Relativistic hydrodynamics itself is an ehective theory that can be used to describe the 

^HTL-resummed calculations of the thermodynamic potential at finite temperature and quark chemical potential(s) 
describe the lattice data well for T > 300 MeV with no free parameters [6-8]. 
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spacetime evolution of the QGP. In the kinetic theory approach to relativistic hydrody¬ 
namics, one obtains the dynamical equations for the bulk variables by taking moments of 
the Boltzmann equation. Ideal hydrodynamics [12-14] and later on viscous hydrodynamics 
[15-46] have been used to study the QGP created in heavy-ion collisions and have proven 
to be quite successful. Recently, anisotropic hydrodynamics [47? -67] has been developed 
in order to extend the range of applicability of relativistic hydrodynamics to situations in 
which the QGP possess a high degree of momentum-space anisotropy (for a recent review, 
see Ref. [68]). 

In most cases, however, the manner in which the EoS is imposed in hydrodynamics is 
somewhat uncontrolled. In many cases, one derives the hydrodynamic equations for a con¬ 
formal system and then imposes an EoS to relate the components of the energy-momentum 
tensor. We refer to this as the “standard EoS” method. However, since QGD is a non- 
conformal theory with a running coupling constant that depends strongly on the temperature 
near Tg, it is more self-consistent to take into account the breaking of conformal invariance 
from the beginning, which results in additional terms in the evolution equations and new 
transport coefficients. Some progress in this direction has been made in the last year, both 
in the context of second-order viscous hydrodynamics [44] and anisotropic hydrodynamics 
[62], however, in both of these previous works, the underlying microscopic picture was that 
of a gas of particles with temperature-independent masses. One would like to incorporate 
the temperature-dependence of the particle masses into the dynamical equations such that 
the equations themselves are consistent with the breaking of conformal invariance and the 
quasiparticle picture at high temperatures. In this paper, we present a method for doing 
this in the context of anisotropic hydrodynamics. Our method is to incorporate the effects 
of a temperature-dependent quasiparticle mass into the Boltzmann equation by taking into 
account extra terms which come from the spacetime gradients of the thermal mass. We 
show that adding the necessary additional term to the Boltzmann equation and enforcing 
energy-momentum conservation require one to introduce a non-equilibrium background field 
Bg^'^ to the energy-momentum tensor as was found by previous authors [10, 11], e.g. 

r"" = + Bj'*" . (1) 

This extra background contribution can be shown to become precisely the additional term 
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necessary to enforce thermodynamic consistency in the equilibrium limit, however, in prac¬ 
tice, we allow it to be a non-equilibrium quantity. 

The new method above will be referred to herein as the “quasiparticle EoS”. We com¬ 
pare results obtained using this method to results obtained using the canonical method for 
imposing a realistic equation of state. For this purpose, we reduce the dynamical equations 
in both cases to those appropriate for 0-l-ld boost-invariant and transversally-homogeneous 
expansion subject to a relaxation-time approximation collisional kernel. With this setup, 
comparisons of the evolution of the effective temperature, pressure anisotropy, and bulk 
correction to the pressure for different values of the shear viscosity to entropy density ratio 
are presented for both isotropic and anisotropic initial conditions. We demonstrate that 
the temperature evolution obtained using the two EoS methods is nearly identical and that 
there are only small differences in the pressure anisotropy. However, we hnd that there are 
signihcant differences in the evolution of the bulk pressure correction, which could poten¬ 
tially be important for determining the correct form of the particle distribution function on 
the freezeout hypersurface in phenomenological applications. 

The structure of the paper is as follows. In Sec. II, we present the notation and conventions 
we use in the paper. In Sec. Ill, we review the necessary setup including the anisotropic dis¬ 
tribution function, basis vectors necessary in different cases, and the lattice-based equation 
of state we will use. In Sec. IV, the Boltzmann equation and its generalization to quasiparti¬ 
cles with temperature-dependent masses is discussed. In Sec. V, we take different moments 
of distribution function in order to derive expressions for the particle current, energy den¬ 
sity, and components of the pressure. In Sec. VI, the 3-|-ld dynamical equations for mas¬ 
sive anisotropic hydrodynamics are derived and then simplihed assuming boost-invariance 
together with either cylindrical-symmetry or transversally-homogeneity. In Sec. VII, we 
obtain the 0-l-ld dynamical equations for the quasiparticle EoS and standard EoS cases. 
In Sec. VIII, our numerical results obtained using both methods for a boost-invariant and 
transversally-homogeneous system are presented. Sec. IX contains our conclusions and an 
outlook for the future. All necessary identities and function dehnitions are collected in 
Apps. A-B. 
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II. CONVENTIONS AND NOTATION 


A parentheses in the indices indicates a symmetrized form, e.g. + A^^)/2. 

The metric is taken to be “mostly minus” such that in Minkowski space with = (t, x, y, z) 
the line element is ds"^ = g^ydx^dx'' = dt^—dx‘^—dy'^—dz‘^. We also make use of the transverse 
projector, — u^u''. When studying relativistic heavy-ion collisions, it is convenient 

to transform to variables dehned by r = y/t^ — z"^, which is the longitudinal proper time, and 
^ = tanh“^(; 2 /t), which is the longitudinal spacetime rapidity. If the system is additionally 
cylindrically symmetric with respect to the beam-line, it is convenient to transform to polar 
coordinates in the transverse plane with r = y/x‘^+y‘^ and 0 = tan~^{y/x). In this case, the 
new set of coordinates x^ = (r, r, 0, dehnes polar Milne coordinates. Finally, the invariant 
phase space integration measure is dehned as 


dP = Ndof 


d^p 1 
{27ifE 



( 2 ) 


where Vdof is the number of degrees of freedom and N = Ndoi/{2n)^. 


III. SETUP 

In this paper, we derive non-conformal anisotropic hydrodynamics equations for a sys¬ 
tem of quasiparticles with a temperature-dependent mass. To accomplish this goal, an 
effective Boltzmann equation for thermal quasiparticles is obtained. We then take moments 
of the resulting kinetic equation to obtain the leading-order 3-|-ld anisotropic hydrodynam¬ 
ics equations. Using a general set of basis vectors, the equations are expanded explicitly 
and then various simplifying assumptions (e.g. boost-invariance, etc.) are imposed to re¬ 
duce the equations from their general 3-|-ld to a 0-l-ld form appropriate for describing a 
boost-invariant and transversally-homogenous QGP. The obtained 0-l-ld equations are then 
solved numerically for our tests, however, the method used to obtain the 3-|-ld leading-order 
anisotropic hydrodynamics equations can be used without lack of generality. 
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A. Basis Vectors 


A general tensor basis can be constructed by introducing four four-vectors which in the 
local rest frame (LRF) are 


'^LRF = (1) 0; 0) 0) ! 

"^LRF = 0, 0) , 

^LRF = 0) , 

^LRF = (0)0)0)1)• (3) 

One can dehne the general basis vectors in the lab frame (LF) by performing the Lorentz 
transformation necessary to go from LRF to the LF. The transformation required can be 
constructed using a longitudinal boost ■d along the beam axis, followed by a rotation tp 
around the beam axis, and hnally a transverse boost by d±_ along the x-axis [52, 53]. This 
parametrization gives 

= (cosh 6± cosh ■d, sinh 6± cos <p, sinh 6j_ sin (p, cosh 6j_ sinh '&) , 

= (sinh d_L cosh cosh 6± cos cosh 6± sin <p, sinh 6± sinh "d), 

= (0, — sin (p, cos p, 0), 

= (sinh "d, 0, 0, cosh'd), (4) 

where the three helds'd, p, and 9± are functions of Cartesian Milne coordinates (r, x,|/, <^). 
Introducing another parametrization by using the temporal and transverse components of 
flow velocity 


Mo = cosh dj_ , 
Ux = u± cos p , 
Uy = u± sin p , 


(5) 

( 6 ) 
(7) 
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where u±_ = ^Jul. + Uy = ^/uq — 1 = sinh6*_L, one has 

= {uq cosh Ux,Uy, Uq sinh ■^), 

= (u± cosh i!), ^ ^ u± sinh li ), 

V u± u± / 

(o,-^,—,o), 

V u± u± / 

= (sinh'i?, 0, 0, cosh'i?). (8) 

For a boost-invariant and cylindrically-symmetric system, one can simplify the basis vectors 
by identifying '& = <:; and (p = 4> where ? and 0 are the spacetime rapidity and the azimuthal 
angle, respectively. In this case, the basis vectors (4) simplify to 

= (cosh 6j_ cosh sinh 6j_ cos 0, sinh 6j_ sin cj), cosh 6j_ sinh <^), 

= (sinh 6± cosh cosh 6± cos 0, cosh 6± sin 0, sinh 6± sinh <^), 

= (0, — sin 0, cos 0, 0), 

= (sinh <^, 0, 0, cosh . (9) 

In the case of a transversally-symmetric system, the transverse flow u± is absent, i.e. 6± = 0, 
and, as a consequence, one has 


= (cosh 0,0, sinh <^), 

= (0, cos 0, sin 0,0), 

= (0, — sin 0, cos 0,0), 

= (sinh 0, 0, cosh <^). (10) 

Note that in the last case, X^ and Y^ are simply unit vectors pointing along the radial and 
azimuthal directions, respectively. 
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B. Ellipsoidal form including bulk pressure degree of freedom 


In the non-conformal case, anisotropic hydrodynamics is dehned throngh the introdnction 
of an anisotropy tensor of the form [53, 62] 

, ( 11 ) 

where is four-velocity, is a symmetric and traceless tensor, and $ is associated with 
the bulk degree of freedom. The quantities u^, and d* are understood to be functions 
of spacetime and obey = 1, = 0, = 3, and = 0; therefore, one has 

= 1 — 3$. At leading order in the anisotropic hydrodynamics expansion one assumes 
that the one-particle distribution function is of the form 

, ( 12 ) 

where A has dimensions of energy and can be identihed with the temperature only in the 
isotropic equilibrium limit = 0 and <I> = 0).^ We note that, in practice, /iso need 
not be a thermal equilibrium distribution. However, unless one expects there to be a non- 
thermal distribution at late times, it is appropriate to take /iso to be a thermal equilibrium 
distribution function of the form fiso{x) = /eq(x) = (e^^ -|- a)“^, where a = ±1 gives Fermi- 
Dirac or Bose-Einstein statistics, respectively, and a = 0 gives Boltzmann statistics. From 
here on, we assume that the distribution is of Boltzmann form, i.e. a = 0. 

C. Dynamical Variables 

Since the most important viscous corrections are to the diagonal components of the 
energy-momentum tensor, to good approximation one can assume that = diag(0,^) with 
^ and = 0. In this case, expanding the argument of the square root appearing 

on the right-hand side of Eq. (12) in the LRF gives 

^Herein we assume that the chemical potential is zero. 


f{x,p) = /is 




(13) 





where i G {x, y, z} and the scale parameters ai are 


a. 


+ + 


( 14 ) 


Note that, for brevity, one can collect the three anisotropy parameters into vector ot = 
{ax,(yy,(yz)- In the isotropic eqnilibrinm limit, where = $ = 0 and ai = 1, one has 
= (p ■ and A —)■ T and, therefore. 



( 15 ) 


Ont of the fonr anisotropy and bnlk parameters there are only three independent ones. In 
practice, we nse three variables as the dynamical anisotropy parameters since, by nsing 
Eq. (14) and the tracelessness of one can write <1> in terms of the anisotropy parameters, 
® In the transversally-homogeneons case, one has a^ = ay and, as a resnlt, 

there are two independent anisotropy parameters. Note that, for conformal systems, one 
has 4) = 0 and in this case there are then only two independent anisotropy parameters in 
3+ld. 

D. Equation of state 

Herein, we consider a system at finite temperatnre and zero chemical potential. At 
asymptotically high temperatnres, the pressnre of a gas of qnarks and glnons approaches 
the Stefan-Boltzmann (SB) limit, Vsb = = Adof7"^/7r^ = tt^T^ — 1 + \NcNf) /45. 

We will take W = A^/ = 3 in what follows. At the temperatnres probed in heavy-ion colli¬ 
sions there are important corrections to the SB limit and at low temperatnres the relevant 
degrees of freedom change from qnarks and glnons to hadrons. The standard way to de¬ 
termine the QGP EoS is to nse non-pertnrbative lattice calcnlations. For this pnrpose, we 
nse an analytic parameterization of lattice data for the QCD interaction measnre (trace 
anomaly), Aq = Teq — 3Peq, taken from the Wnppertal-Bndapest collaboration [69] 


/eq(T) 


1 -F g-^t + g2t‘^ 


h _^ /o[tanh(/it-F /s) l] 



(16) 
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(b) 




FIG. 1. Panel (a) shows the energy density and pressure scaled by their respective Stefan- 
Boltzmann limits as a function of temperature. Panel (b) shows the speed of sound squared 
as a function of temperature. 


with t = T/ (0.2 GeV). For Nf = 2 + 1 (2 light quarks and one heavy qnark) the parameters 
are = 0.1396, hi = -0.1800, hs = 0.0350, /o = 2.76, /i = 6.79, /s = -5.29, gi = -0.47, 
g 2 = 1.04, and = 0.01. 

The pressnre can be obtained from an integral of the interaction measnre 


P=,(r) 

rp/l 


dTh^{T) 

rp 


( 17 ) 


where we have assnmed Veq{T = 0) = 0. Having Peq(T), one can obtain the energy density 
£’eq using Seq(T) = 3Veq(T) + Jeq(T). In the limit T —)■ cx), the system tends to the ideal 
limit as expected.^ The temperature dependence of the resulting equilibrium energy density, 
pressure, and speed of sound squared (c^ = dVeq/dSeq) are shown in the two panels of Fig. 1. 


Method 1: Standard equation of state 


In the standard approach for imposing a realistic EoS in anisotropic hydrodynamics, one 
derives the necessary equations in the conformal limit and exploits the conformal multi¬ 
plicative factorization of the components of the energy-momentum tensor [47, 48]. With 
this method, one relies on the assumption of factorization even in the non-conformal (mas¬ 
sive) case. Such an approach is justified by the smallness of the corrections to factorization 


^In the original parametrization presented in Ref. [69] the authors used /13 = 0, however, as pointed out in Ref. [ 66 ], 
choosing /13 = 0 gives the wrong high temperature limit. 
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in the massive case in the near-equilibrium limit [66]. For details concerning this method, 
we refer the reader to Refs. [54, 66]. Although this method is relatively straightforward 
to implement, it is only approximate since for non-conformal systems there is no longer 
exact multiplicative factorization of the components of the energy-momentum tensor. This 
introduces a theoretical uncertainty which is difficult to quantitatively estimate. 

Method 2: Quasiparticle equation of state 

Since the standard method is only approximate, one would like to hnd an alternative 
method for imposing a realistic equation of state in an anisotropic system that can be ap¬ 
plied for non-conformal systems. In order to accomplish this goal, we implement the realistic 
EoS detailed above by assuming that the QGP can be described as an ensemble of massive 
quasiparticles with temperature-dependent masses. As is well-known from the literature 
[9], one cannot simply substitute temperature-dependent masses into the thermodynamic 
functions obtained with constant masses because this would violate thermodynamic consis¬ 
tency. For an equilibrium system, one can ensure thermodynamic consistency by adding a 
background contribution to the energy-momentum tensor, i.e. 

rC = re,.,,,,, + , ( 18 ) 

with Beq = Beq{T) being the additional background contribution. The kinetic contribution 
to the energy momentum tensor is given by 

TZetic,eq = j dP p^p’^ feqZ, p) . (19) 

For an equilibrium Boltzmann gas, the number and entropy densities are unchanged, 
while, due to the additional background contribution, the energy density and pressure are 
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(a) 


(b) 




FIG. 2. In panel (a) we plot the temperature dependence of the quasiparticle mass scaled by the 
temperature obtained using Eq. (26). In panel (b) we plot the temperature dependence of the 
background term B^q scaled by the temperature obtained using (25). 

shifted by +i?eq and —B^q, respectively, giving 


neq(T,m) 

= AtiNT^ ^eq-^2 (^ 

eq) 5 

Seq{T,m) 

= AuNT^ml^ 

AK 2 ( 

rheq) + 

£eq{T,m) 

= AirNT^ml^ 

eq 

3K2 

rheq) + 

VeqiT^m) 

= ATTNT^ml^K 2 {m 

eq) B 


"eq-* 


"eq-* 


"eqy 


"eqy 


+ Beq , 


( 20 ) 

( 21 ) 

( 22 ) 

( 23 ) 


where m^q = m/T with m implicitly depending on the temperature from here on. In order 
to fix Beq, one can require, for example, the thermodynamic identity 


TSeq — Sea + ^ea — T 


dR 


eq 


-eq 


eq 


dT ’ 


( 24 ) 


be satisfied. Using Eqs. (22), (23), and (24) one obtains 


dB, 


eq 


1 


dT 


2 dT 


dP feq{x,p) 
~ n , , dm 


(25) 


If the temperature dependence of m is known, then Eq. (25) can be used to determine Beq. 
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In order to determine m, one can use the thermodynamic identity 


+ ^eq = T^eq = dvriVT'^ (rheq) . (26) 

Using the lattice parameterization (16) to compute the equilibrium energy density and pres¬ 
sure, one can numerically solve for m{T). In Fig. 2a, we plot the resulting solution for 
m/T as a function of the temperature. Once m is determined using Eq. (26), one can solve 
Eq. (25) subject to the boundary condition Beq{T = 0) = 0 to hnd i?eq(F). We note that, 
using this method, one can exactly reproduce the lattice results for energy density, pressure, 
and entropy density. In Fig. 2b, we plot the resulting solution for the normalized quantity 
Beq{T)/T‘^ as a function of the temperature. 


IV. BOLTZMANN EQUATION AND ITS MOMENTS 

In this paper, we derive the necessary hydrodynamical equations by taking the moments 
of Boltzmann equation. In what follows, we specialize to the case that the collisional kernel is 
given by the relaxation-time approximation (RTA), however, the general methods presented 
here can be applied to any collisional kernel. If the particles that comprise the system have 
temperature-independent masses then the Boltzmann equation is of the form 

Jfdj = -Clf]. (27) 

The function C[f] at right-hand side of the equation is the collisional kernel containing all 
interactions involved in the dynamics. In RTA, one has 

C[/] = ^(/-/eq). (28) 

In this relation, /eq denotes the equilibrium one-particle distribution function (15) and Teq is 
the relaxation time which can depend on spacetime but which we assume to be momentum- 
independent. To obtain a realistic model for Teq, which is valid for massive systems, one 
can relate Teq to the shear viscosity to entropy density ratio. For a massive system, one has 
[70, 71] 

req{T)Veq{T) . 

v{T) = — -- K(meq). (29) 
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In this formula the function k{x) is defined as 


3 K^i^x) 1 Ki{x) 71 1 — xKq{x)L_i{x) — xKi{x)Lq{x) 
'^K2{x) ~x^ K2{x) ~ 2 K2{x) 


where Kn{x) are modified Bessel functions of second kind and Ln{x) are modified Struve 
functions. Assuming that the ratio of the shear viscosity to entropy density, rj/Seq = f], is 
held fixed during the evolution and using the thermodynamic relation + 'Peq = TSeq one 
obtains 


re,{T) = 


ISr^ 


1 + 


£^(T) 


V Peq(r) 


Note that, in the massless limit, m —t 0, one has K(meq) —t 12, giving 


(31) 


re,(T) 


5rj 

4Pe,(T) ■ 


(32) 


A. Effective Boltzmann Equation 

If the quasiparticles have a temperature-dependent mass, one has to generalize the Boltz¬ 
mann equation in order to take into account gradients in the mass. Generally, the Boltzmann 
equation for on-shell quasiparticles can be written as [10] 

fd dE d dE , rdf\ 

\dt^ dp (9x (9x dp) V^^/coii 

where the “external force term” —dE/d'K does not vanish in the thermal mass case since the 
particle energy depends on the mass and hence on the local temperature of the system. As 
a result, the temperature-dependence of the mass acts as an external force in the dynamics. 
Using the on-shell energy relation E = ^Jp^^^rw? and defining the collisional kernel as 

(34) 

in covariant form one has 

P^d^f + ]^dim^d\p)f = -C[f] , ( 35 ) 

where = ( \Jp‘^ + m?, p) is the on-shell momentum four-vector, z is a spatial coordinate 
index, and = —d/dp^. Note that the extra term, {dim?/2)d]^^-^f, corresponds precisely 
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to the result obtained from derivation of the Boltzmann equation using quantum held the¬ 
oretical methods [72], 

We mention that an alternative way of deriving the effective Boltzmann equation above 
can be found in a recent paper of Romatschke [11]. In this paper, a general Boltzmann 
equation for off-shell particles (with constant mass) is hrst derived using the evolution of 
a single particle distribution function along eight-dimensional phase space geodesics, where 
the possibility of curved spaces is taken into account using geometrical covariant derivatives. 
Then, by adding a temperature-dependent background term to temperature-dependent 
masses are taken into account in a way that guarantees both thermodynamic consistency in 
the equilibrium limit and energy-momentum conservation, in general. Finally, the on-shell 
version of the effective Boltzmann equation for quasiparticles with a temperature-dependent 
mass in a general curved space time is obtained. The hat spacetime limit of the ehective 
Boltzmann equation derived by Romatschke is the same as Eq. (35). 

B. Moments of Boltzmann Equation 

If one is interested in the evolution of the bulk properties of a system, one can use low- 
order moments of the Boltzmann equation. By calculating moments of Boltzmann equation 
one obtains evolution equations for tensors of diherent ranks, with the hrst moment giving an 
evolution equation for the energy-momentum tensor and the second-moment describing the 
evolution of a rank three tensor. Taking the zeroth, hrst, and second moments of Boltzmann 
equation gives, respectively 


= - 

[dPC[f], 

(36) 

= - 

j dPp''C[f \. 

(37) 

= - 

1 dPpy^cif] , 

(38) 


15 


where the particle four-current energy-momentum tensor and the rank-three tensor 
-£^iu\ given by 


= 

^ dPp^f{x,p ), 

(39) 

rjnfil/ _ 

l‘dPp>^p'^fix,p) + Bg^^, 

(40) 


jdP p^p'^p^f{x,p). 

(41) 


We note that we have introduced the non-equilibrium background held B = B{cx., A), which 
is the analogue of the equilibrium background B^q in order to guarantee that the correct 
equilibrium limit of is obtained. In the process of the derivation one hnds that, in order 
to write the energy momentum conservation in the form given in Eq. (37), there must be a 
differential equation relating B and the thermal mass 

d^B = - jdPf{x,p). (42) 

In practice, one can use (42) to write the derivative of B with respect to any variable in 
terms of the derivative of the thermal mass times the E~^ moment of the non-equilibrium 
distribution function. 

V. BULK VARIABLES 

In this section, bulk variables, i.e. number density, energy density, and the pressures, are 
calculated by taking the projections of and 

A. Particle current 4-vector 

The particle current four-vector = (u, J) is dehned in Eq. (39). One can expand 
using the basis vectors as 


. 


(43) 
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Using Eqs. (13) and (39) one has 


= (n, 0) = nu^ , (44) 

where n = aneq(A,m) and a = axQya^- 

B. Energy-Momentum Tensor 

The energy-momentum tensor is dehned in Eq. (40). Expanding it using the basis 
vectors one obtains 


+ VxX^X^ + VyY^^Y’' + . (45) 

Using Eqs. (13), (40), and (45) and taking projections of one can obtain the energy 
density and the components of pressure 

E = 7 ^ 3 ( 0 :, m) + i? , 

Vx = 'HsxioL^m) A^ - S , 

Vy = T-isyict, - B , 

V,='H3L{c^,m)X^-B, (46) 

where m = m/X. In the transversally-symmetric case one has Vt = Vx = Vy and Vl = Vz 
and Eq. (45) simplihes to 


= {E + Vt) - Vrg^'' + (Vl - Vt) Z^^Z’' , (47) 


where 

E = 773 ( 0 :, 7 %) X‘^ + B , 

Vt = 'f(-‘iT(oLi fXi) A^ — B , 

VL='H^L(cc,m)X^-B. (48) 

The various "H-functions appearing above are dehned in App. B 1. 
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VI. DYNAMICAL EQUATIONS 


In order to obtain the dynamical equations from Eqs. (36)-(38), one needs the tensor de¬ 
composition of and using the basis vectors. Herein the general 3-|-ld equations 

for a system with temperature-dependent masses are obtained in the RTA. We then simplify 
to the case of 0-l-ld transversally-symmetric case by the taking necessary limits. In what 
follows, the convective derivatives and divergences 9a, with a G {u,x,y, z}, are dehned 
in App. A. 


A. Zeroth moment 


The evolution equation for the particle four-current (36) in the RTA is 


= — (Ueq - n). 
Aq 


Using Eq. (44) one has 

Dy^n + nOy = — (ueq - n). 

^eq 

In the case of 0+ld, this simplifies to 


drU -I-= -(Ueq — u) . 

'T Aq 


(49) 


(50) 


(51) 


B. First Moment 

The conservation of energy and momentum is enforced by = 0. This requires that 

both the left and right hand sides of Eq. (37) vanish. The vanishing of the right-hand side 
of this equation results in a constraint equation that can be used to write T in terms of the 
non-equilibrium microscopic parameters a and A. Using (13), (15), and (28) one obtains 

4inetic = kinetic,eq, or more explicitly 


nj \4 _ nj rpA 

ttsA — rt^^eq-L ■ 


(52) 


Turning to the left hand side, using Eq. (45) and taking U-, X-, Y-, and X-projections, 
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one obtains four independent equations 


Du£ + £-Su + 'PxUfxDxX^ + T’yiinDyY^ + — 0 , 

DxV:, + - SXyD^U^ - VyXyDyY^ - V,XyD,Z>^ = 0 , 

DyVy + VyOy - SY^D^u^ - VxY^D^X^ - P,YyD,Z>^ = 0 , 

DzVz + - SZ^DuU^ - V^Z^D^X^^ - VyZyDyY>^ = 0 . (53) 

In the 0+ld case, using Vt = Vx = Py and Pl = Pz and taking the appropriate limits, as 
explained in App. A, one can simplify Eqs. (53) to 

^ £ + Pl 

dr£ = -, (54) 

T 

drPT = O^Pt = d,PL = 0. (55) 

Eqs. (55) are consequences of boost invariance and transverse homogeneity in the 0+ld case 
and, as a result, the only independent dynamical equation is Eq. (54). 

C. Second moment 

The second moment of Boltzmann equation in the RTA is 

_ j{uqX)^2 ^ }_ ^ 

where can be obtained from Eq. (41) by taking / —)• /eq. For a distribution function 
of the form specihed in Eq. (13), the only non-vanishing terms in Eq. (41) are those with 
an even number of similar spatial index. As a result, one can expand over the basis 
vectors as 


X = X„ [m (g) M (g) u] 

+ X3,[MgXgX + XgMgX + XgXgM] 

+ Xy[ugFgF + ygugF + ygygu] 

+ Xz[u ® Z ® Z + Z ® u ® Z + Z ® Z ® u]. (57) 
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Evaluating the necessary integrals using the distribution function (13), one hnds 

Xu= ( y^^a^)aZeci{X,m) + am^raeq(A,m), (58) 

i 

Xi = a Xeq(A, m) , (59) 


where 

Xcq(A,m) = AttNX^771^X3(7%). 

Note that, in general, one has Xu—Yhi ~ am^neq(A, m) and limm-s>o Xu 
Eq. ( 56 ) and taking its uu-, XX-, YY-, and XX-projections gives 

DuXu + XuOu + XLxU^DxX^ + 2Xy7j,^DyY^ + 2X^77 ^DzZ^ — 72 Du7ti^ = 

DuXx + XxiOu + 2uyDxX^) = 
DuXy + Xy(9u + 2,7J,yDyY^) = 
DuXz + Xz{0u + 2uyDzZ^) = 

Also, taking uX-, uE-, and uX-projections one can End 

DxXx + Xx9x + (Xx + Xu)UyDuX^ — XyXyDyY^ — XzXyDzZ^ — 
DyXy + XyOy X {Xy X Xu)U yD uY ^ - X.E^D.X^ - XzYyD z Z ^ - 
X>zXz + Xz6z + (X^ + Xu)UyDuZ^ — XxZyDxX^ — XyZyDyY^ — 

and hnally projecting with XE, XX, and EX gives 

Xx{YyDuX>^ + YyDxU^) + Xy{XyDuY^ + XyDyU^^) = 
Xx{ZyDuX^ + ZyDxU>^) + XziXyDuZf^ + XyDzU^) = 
Xy{ZyDuY^ + ZyDyU^^) E X,(E^X„X'^ E YyDzU^^) = 


(60) 

= Y.iXi- Expanding 

— (X„,eq - Xu ), (61) 
Eq 

— (Xeq-Xj, (62) 

— (Xeq-Xj/), (63) 

Eq 

— (Xeq-X,). (64) 

^eq 

^TlDxTTl'^ = 0 , (65) 

^nDy'm'^ = 0 , (66) 

^nDzTJi'^ = 0 , (67) 

0 , ( 68 ) 
0, (69) 

0. (70) 


It can be shown that Eq. (61) is not independent. One can subtract the sum of Eqs. (62)- 
(64) from it to obtain 

9 

Tfl 

7n^{Du7i + n9u) = —(ueq - n). (71) 

Eq 
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This equation is the same as Eq. (50) for non-vanishing mass. In the 0-l-ld case, one has 
Xx = Xy and Eqs. (62)-(64) simplify to 

drlogI,. + - = X{^.l\ (72) 

8,logI. + - = T('52_i'), (73) 

^eq X^ / 

Finally, we note that Eqs. (65)-(70) are trivially satished in the 0-l-ld case. 

D. Selection of relevant eqnations of motion 

For the 0-l-ld case, we need four equations for the four independent parameters. A, T, a^- 
Using the equations derived thus far up to the second moment of Boltzmann equation, we 
have five independent equations. Herein, we use the equations obtained solely from the first 
and second moments of the Boltzmann equation which give Eqs. (52), (54), (72), and (73).^ 

VII. 0-hlD DYNAMICAL EQUATIONS 

In this section, we present the dynamical equations for the “quasiparticle EoS” and the 
“standard EoS” cases. For simplicity, we present only the 0-l-ld case herein. We postpone 
the 3-|-ld numerical comparisons to a future work. 

A. Quasiparticle equation of state 

One potential complication encountered when using temperature-dependent masses is 
that the first moment equation will involve the background contribution B and its proper¬ 
time derivative, since on the left-hand side of (54) one has the total energy density which 
includes the background contribution. In practice, however, all derivatives of B can be 
written in terms of derivatives of m using Eq. (42). For the 0-l-ld case, we only need the 
proper-time derivative of B. Taking the distribution function to be of the form (13) and 
using Eq. (42) one obtains 

drB = m) drrn^ . (74) 

^For the “quasiparticle EoS” case one obtains quite similar results if one instead uses the equation obtained from the 
zeroth-moment (51); however, in the “standard EoS” case, one finds that using the zeroth moment equation (51) 
results in solutions that do not approach the isotropic equilibrium limit at late times. 
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In this way, all proper-time derivatives of B necessary for the evolntion eqnations can be 
obtained from derivatives of the thermal mass and knowledge of the non-equilibrium mi¬ 
croscopic parameters which enter the "Has function. However, in order to obtain the total 
energy density or pressures one needs to know B itself. Our procedure will be to integrate 
the dynamical equations to a very late proper time Tf when the system is close to equilib¬ 
rium and then integrate Eq. (74) backwards in time from r/ to the intial time tq subject to 
the boundary condition that B{Tf) = Heq(T(r/)). 

Using Eqs. (48), (59), and (60) one can expand (54), (72), and (73) to obtain 


d B ^ 

Al-i^dr log A Vtmdr log rh Vtidr log log -f H-^ = 0 , 

A r 

Adr log ax + dr log az + bdr log \ + dr log (mbit's (m)) H— 


1 

1 


A773(meq) ■ 

"^eq 

Q 

Q 

_1 

Ia/ 

K^irh) 


3 

2dr log ax + “idr log az + 59^ log \ + dr \og{m^K'i{m)) -\— 


1 

\ 1 

fT^ 

s^K^{rheq) ■ 

'^eq 

alal 

a/ 

Ks{m) 


(75) 


(76) 


(77) 


where ^l-, and are defined in App. B 1. 

One can perform some algebra to change the matching condition (52) into a differential 
equation which is more convenient to solve since we then only have to solve a system of 
coupled ordinary differential equations. Taking a derivative of Eq. (52) with respect to r 
and using Eq. (54), one obtains 


4773,eq^r log T + Om.eq^r log ^eq + ^ = 0 . 


r 


J'A 


(78) 


In all equations above one can use Eq. (31) for Teq. 


B. Standard equation of state 

We now present the details of our implementation of the “standard EoS” method. In this 
case, one takes the particles to be massless, m —)■ 0, and hence H —)■ 0. For the massless 
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transversally-symmetric case, Eqs. (48) become 


£ = Moi)X\ 

Vt = 1 

VL = n^L{oc)\\ (79) 


where all 'H-functions are defined in App. B 2. As we can see from the above eqnations, there 
is a mnltiplicative factorization of the energy density and pressnres into a fnnction that only 
depends on the anisotropy parameters and a fnnction that only depends on the scale A. 
For a massless conformal Boltzmann gas, one has £eq{T) = 247rA^T^ and Veq{T) = SttA^T^. 
Using these relations, one can rewrite Eqs. (79) in terms of the equilibrium thermodynamic 
functions 


8 = 
'Pt — 
Vl = 


i’eq(A) 




3Peq(A) 


4 


a. 


3Peq(A) 


OIxT^2L - 

ax 


(80) 


These formulas suggest that, in order to impose a realistic EoS, one only has to replace 
Teq(A) and Peq(A) by the results obtained from lattice QCD calculations. 

In order to obtain the necessary dynamical equations, one has to take the limit m —)■ 0 
of the equations obtained from the moments of the Boltzmann equation and substitute £ 
and Vt,l from Eq. (80). For the first moment equation, starting from Eq. (54) and using 
Eq. (80) one obtains 


a^log£:eq(A) + (1 + x)9^1oga^ + (3 - x)dAogax, = - - - , (81) 

'^£eq ( A ) 

with X = 71211^2 ■ Taking the limit m —)■ 0 and i? —)■ 0 of the second-moment equations 
(76) and (77), one obtains 


Adr log ax + dr log ar + ^dr log A - 

T 

3 

2dr log ax + 3dr log a^ + 59,- log A -I— 

r 



5 1 

a\az 
5 1 



(82) 

(83) 
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For the matching relation which gives T in terms of the microscopic parameters, one can 
use £^(A) = EecSJ') and Eq. (81) to hnd 


dA0g£e^{T) 


1 3Peq(A) 
r r£:eq(A)^ 


In all equations above one can use Eq. (32) for Teq. 


(84) 


Vm. RESULTS 

In this section, we present the results of numerically integrating the dynamical equations 
using the “standard EoS” and the “quasiparticle EoS” methods. In both cases, we specialize 
to the 0+ld case. We take the initial proper time to be tq = 0.25 fm/c and the final time to 
be Tf = 500 fm/c. In all cases, the initial temperature is taken to be Tq = 600 MeV which 
is appropriate for LHC heavy-ion collisions at ^snn = 2.76 TeV. The hnal time used here is 
very long compared to the timescales relevant for heavy-ion collisions, but we are interested 
in the late-time approach to isotropic thermal equilibrium in both approaches. Additionally, 
as mentioned previously, in order to determine B{t), we solve the the differential equation 
(74) by evolving it backwards in proper time subject to a boundary condition that H(tj) = 
Heq(T(r/)) and, consequently, we should evolve the system to a late proper-time at which 
the system is close to isotropic thermal equilibrium. 

Before proceedings to our results, we need to dehne one quantity which has yet to be 
dehned, namely the bulk correction to the pressure. In viscous hydrodynamics, the energy- 
momentum tensor is expressed generally as 

- (Peq + n) , ( 85 ) 

where £eq = Teq(T) and "Peq = T’eq(T) are the equilibrium energy density and pressure 
evaluated at the effective temperature. In the dehnition above, is the shear tensor 
and n is the (isotropic) bulk correction. Since is a traceless tensor, tt)) = 0, which is 
transverse to the fluid four-velocity, = 0, one hnds that the bulk correction can be 

computed from 

n = - Peq =\i^L + 2Vt) - Peq • (86) 
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pj'Pt T/r( 


(a) (b) 



0.5 1 5 10 50 100 500 0.5 1 5 10 50 100 500 

T [fm/c] r [fm/c] 

FIG. 3. The four panels show: (a) the effective temperature scaled by Tq, (b) times the 

energy density scaled by the initial energy density, Sq, ( c ) the pressure anisotropy, and (d) the bulk 
correction to the pressure scaled by Veq- For this hgure we took 47rr//s = 1. 

For the case of a temperature-dependent mass, one can use Eqs. (23) and (48). For the 
massless case, one can use Eqs. (17) and (80). 


Numerical results 

We now turn to our numerical results. In all plots, we compare the two methods for 
implementing the EoS in anisotropic hydrodynamics. For the curves labeled “Quasiparticle 
EoS”, we solve the dynamical equations specihed in Sec. VIIA and for the “Standard EoS”, 
we solve those in Sec. VIIB. For purposes of the comparison, we match physical quantities 
rather than the microscopic parameters at Tq. In practice, this means that we specify an 
initial temperature Tq, an initial momentum-space anisotropy quantihed by VlaIVta, and 
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PlI'Pt T/r( 


(a) (b) 



FIG. 4. Same as Fig. 3 except here we take A-tit]! s = 3. 


an initial bnlk correction, Ho. In all results figures, we present four panels which correspond 
to: (a) the effective temperature scaled by Tq, (b) times the energy density scaled 

by Eq, (c) the LRF pressure anisotropy, and (d) the bulk correction scaled by the equilibrium 
pressure, Il/Peq- 

In Figs. 3 - 5 we present our results for the case of isotropic initial conditions. In all 
panels, the microscopic parameters were adjusted to achieve = 1 and Hq = 0. 

From panel (a) of this set of hgures, we see that there is excellent agreement between 
the effective temperature predicted by each method for implementing the EoS. In practice, 
we found that, for all initial conditions we considered, the maximum difference between 
the effective temperature obtained using the two approaches was less than on the order of 
1%. To further explore the differences in the “hrst order” quantities, in panel (b) we have 
multiplied the scaled energy density by a factor of (r/ro)"^'^^. If the system behaved as an 
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PlI'Pt T/r( 



FIG. 5. Same as Fig. 3 except here we take A-nr]/s = 10. 


ideal gas undergoing boost-invariant expansion in ideal hydrodynamics, then at late times 
this quantity should approach unity. Any late-time deviations from unity are indicative of 
the corrections to ideal Bjorken scaling. As we can see from panel (b) of Figs. 3-5, the 
energy density evolution obtained using the two approaches is quite close, with the largest 
difference between the two approaches being approximately 4%. 

Considering panel (c) of Figs. 3 - 5, we see that there are larger differences in the pressure 
anisotropy predicted by the two approaches. For this quantity, we see differences as large 
as 20%, however, the behavior of the pressure anisotropy is qualitatively the same overall. 
Finally, we turn to panel (d) of Figs. 3-5 which shows the bulk correction scaled by 
the equilibrium pressure. As we see from these panels, there is a qualitative difference in 
the temporal evolution of the bulk correction when comparing the two approaches. At late 
times, however, both approaches seem to converge to the same asymptotic limiting behavior. 
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PlI'Pt T/r( 


(a) 


(b) 



FIG. 6. Same as Fig. 4 except here we take anisotropic initial conditions. 


Note that the differences in the pressnre anisotropy and bulk correction are already self- 
consistently taken into account in the evolution of the temperature/energy density. In this 
sense, despite having differences in the viscous corrections, the hrst order quantities seem 
to be quite insensitive to whether one uses the quasiparticle EoS method or the standard 
EoS method. That being said, the differences seen in panels (c) and (d) could manifest 
themselves as differences in the particle spectra computed along the hypersurface if these 
two methods are applied to QGP phenomenology. 

Finally, in Fig. 6 we present the same four panels, but in the case of an anisotropic initial 
condition with an oblate momentum-space anisotropy. As can be seen from Fig. 6, even for 
anisotropic initial conditions, the two EoS methods agree extremely well for the evolution 
the effective temperature and energy density. However, similar to the case of isotropic initial 
conditions, we see somewhat larger differences in the evolution of the pressure anisotropy and 
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qualitative differences in the evolution of the bulk pressure correction. For both quantities, 
we see that the two methods have the same late-time asymptotic behavior. Finally, we note 
that the behavior seen in Fig. 6 is indicative of the results we obtained for a variety of 
different non-equilibrium initial conditions. 

IX. CONCLUSIONS AND OUTLOOK 

In this paper, we presented a new method for imposing a realistic EoS in the context 
of anisotropic hydrodynamics. The method relies on a quasiparticle picture of the QGP, 
which is conceptually consistent with the kinetic theory method used to derive the re¬ 
quired hydrodynamic evolution equations from the Boltzmann equation. We discussed the 
fact that the introduction of a temperature-dependent quasiparticle mass requires an addi¬ 
tional background contribution to the energy-momentum tensor. We showed that requiring 
energy-momentum conservation results in a constraint equation on the background contribu¬ 
tion which reduces to the constraint necessary to enforce thermodynamic consistency in the 
isotropic equilibrium limit as found by previous authors [10, 11]. When solving the result¬ 
ing dynamical equations, we allowed the background contribution to be a non-equilibrium 
quantity. This was necessary to self-consistently implement the constraint equation. 

By numerically solving the resulting dynamical equations in the 0-l-ld, we compared 
the results obtained using the quasiparticle EoS method with those obtained using the 
standard method for imposing a realistic EoS in anisotropic hydrodynamics. We found that 
the temperature evolution obtained using the two methods was nearly identical and that 
there were only small differences in the pressure anisotropy. However, we found that there 
were large qualitative differences in the evolution of the bulk pressure correction. These 
conclusions were supported by the presentation of results for both isotropic and anisotropic 
initial conditions and also for different values of the shear viscosity to entropy density ratio, 
however, we internally checked a much larger set of initial conditions/parameter sets and 
found that these conclusions were generic. We note, however, that the difference in the 
bulk pressure correction found does not necessarily imply large corrections for heavy-ion 
phenomenology. As we have shown, first order quantities like the energy density are the 
same to within a few percent when comparing the two approaches. That being said, the 
differences in the bulk pressure in particular could be important when hxing the form of the 
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distribution function on the freezeout hypersurface. 

Looking forward, in a future work we will present numerical comparisons of the two 
approaches beyond the simple case of 0+ld expansion considered herein. Additionally, it 
would be quite interesting to apply the quasiparticle EoS method to obtain the dynamical 
evolution for a non-conformal system with temperature-dependent masses within the context 
of second-order viscous hydrodynamics. Finally, we note that it may be possible to construct 
exact solutions of the RTA Boltzmann equation for a system of particles with temperature- 
dependent masses using methods similar to those in Refs. [73-75]. Additionally, for the case 
of quasiparticle masses that are linear in the temperature, it may be possible to exactly 
solve the RTA Boltzmann equation subject to Gubser flow similar to Refs. [76, 77]. 
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Appendix A: Explicit formulas for derivatives 

In this section, hrst we introduce the notations used in derivation of the general moment- 
based hydrodynamics equations and then, by taking the appropriate limits, we simplify them 
for the transversally-homogeneous 0-l-ld case. Using the definitions 


V = cosh('d — <;)dr + - sinh('d — <i)d ^, 

T 

V = sinh('d — -f - cosh('d — c,)d ^, 

r 

V_L ■ u_L = + dyUy , 

ll_i_ • V j_ = Ux^x T : 

ll_L X V _1_ = Xlxdy Uydx , 


(Al) 


(A2) 
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and four-vectors defined in Eq. (8) one has 


Du = u^df, = UqV + uj_-V±, 

D, = X>^d, = u^V + ^{u^-V^), 

u± 

Dy = Y^^d, = —{u^xV^), 
u± 

D, = Z>^d^ = V . (A3) 


The divergences are defined as 

6u = = Vuq + uo'b'd + V_L ■ u_L, 

e, = = Vu^ + uYD'd + —(Vx ■ ux) - ^(ux ■ Vx)m± , 

M_L UoU\ 

By = dyY^ = —^(u_L ■ v±)(^, 

U± 

9, = dyZf^ = V^, (A4) 

where ip = ta.n~^{uy/ux). 

D(^u± , 

Mo 

UyDoY^ = u±DaP >, 

UyDaZ^ = UoDa'B , 

XyDoX^ = UoDaP^ , 

XyDaZ^ = U±Da^ , 

YyD^Z^ = 0, (A5) 

where a G {u, x, y, z}. Note that contractions such as X^D^Uy are also non-vanishing, how¬ 
ever, such terms can be written in terms of the expressions above by using the orthogonality 
of the basis vectors, i.e. Da{X^Uy) = 0 implies that X^D^Uy = —UyDaX^. 
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1. Simplification for 1+ld 


In the case of boost-invariant and cylindrically-symmetric flow one has ^ cj) and "d — 
where is the spatial rapidity. Using u± = sinh^^, one can rewrite (A 2 ) as 


V = dr, 
V=-d,, 

T 

V7 O , “-L 

V_L ■ U_L = OrU± H-, 

r 

U_L ■ V_L = U±dr , 

rv “-L Q 
U_L X V_L = —. 


Also, the identities in (A 4 ) become 


(A 6 ) 


Du = cosh 9 ±dr + sinh 9 ±dr , (A 7 ) 

= sinh 9 ±dr + cosh 9 ±dr , (A 8 ) 

Dy = ^ 84 ,, (A 9 ) 

D, = -d,, (AlO) 

T 

9 u = cosh 6 *_L -I- dr 9 ±^ + sinh^^ -|- dr 9 ±^ , (All) 

9^0 = sinh 9 ± -|- dr 9 ±^ + cosh 9 ±(^- + dr 9 ±^, (A 12 ) 

0 ^ = g^ = 0 _ (A 13 ) 

In this limit, the only non-vanishing terms in (A 5 ) are 


UyDuX^ = Du9± , 

= D^9x ., 

UyDyY^ = - sinh 9± , 

UnDzZ^ = - cosh6*_L, 

T 

XyDyY^^ = ^cosh0^, 

XyD.Z^^ = i sinh 9^_. (A14) 
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2. Simplification for 0+ld 


For this case, one has 6*^ = 0 and 


In this limit, the only non-vanishing 


Du = dr, 


Dx = dr, 

D 

r 

T 

0u = - , 
T 

9x = - , 


r 

9y = 9, = 0. 

terms in (A5) are 


= - , 
T 

= I ■ 


(A16) 

(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 


Appendix B: special functions 


In this section, we provide dehnitions of the special functions appearing in the body of 
the text. We start by introducing 

/•i 


K2{y, z)=y dx y/ {y‘^ — l)x‘^ + z'^ + 1 


'-1 

y 




+ 1) tanh 


-1 


y^-i 

y2 _|_ ^2 


+ Viy'^ - 1 )( 2 /^ + 2 ;^) 


(Bl) 


'H2T{y,z) = y 


dx(l — 


-1 


y/ — l)x‘^ + + 1 


y 


(2/2 - 1)3/2 


{z^ -I- 2y^ — l) tanh ^ 


y^-i 

y2 _|_ ^2 




(B2) 
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U2L{y,z) = 


dxx^ 


-1 


a/ {y‘^ — l)x‘^ + 2:^ + 1 


A /17/2 1 ir7/2 -J- r2i 

(M 1 lUanh"^ i ^ 

V IP ■‘-IIP ^ ) 

\ -L J Lclilll A / 0 0 

V p + 


(^2 _ 1)3/2 

Derivatives of these functions satisfy the following relations 


(B3) 


dn2{y,z) 

dy 

d'H2{y,z) 

dz 


1 

y 

1 

2 : 


'H2{y,z) +'H2L{y,z)^, 

'H2{y, z) - 'H2L{y, z) - 'H2T{y, z) 


(B4) 

(B5) 


1. Massive Case 


The "H-functions appearing in the dehnitions of components of the energy-momentum 
tensor are 


Uz{cx,m) 

n3x{a,m) 

n3y{a,m) 

n 3 T{a,m) 

'H3L{cx,m) 

n3m{o^,m) 

n3B{a,m) 

f2T(CK, hi) 

Qrn{a,m) 


r-27r 


N OixOiy 


dcj) ai / dpp^feq (yp'^ + 7/2 f 


ttz m 


a± a±p 


Malay J dp cos^ p J dpp^feqiyp"^ + '^2 t ^^, 

Naxal [ dpsin^p j dp p^ feqiypP^r^'^'H2T , 


m 


\q;_l q;_lP 
m 


'H3x{a,m) -h 'H3y{a.,m 


r-27r 


Naxay 


dp a'i J dpp^ /eq (yp^ + m"^^ 'H 2 L ^ ^, 
/eqfVP^ +hl 2 


\q;_l q;_lP 


a^ m 


NaxayWp / dpal_ dpp^ 


a/p^ + 


a± a±p 


n 2 {—,—. 

\q;j_ q;j_p 


^27r poo 


dp / dp p/eq (yp‘^ + m^^ 'H 2 B f ^, 


Maxay I 

Jo JO 

JJ-3 + JJ-OT , 

JJ-0 + JJ-OL ) 

JJ-O — JJ-OL — ‘JJJ-3T — JJ-3m , 


a^ m 


q;j_ q;j_p 


(B6) 

(B7) 

(B8) 

(B9) 

(BIO) 

(Bll) 

(B12) 

(B13) 

(B14) 

(B15) 


34 



where aj_ = cos^ (p + ay sin^ p and 


'H 2 B{y, z) = U 2 T{y, z) + ^ tanh ^ W ^ ^ . (B16) 

y^ Vi/2 - 1 y^ + 

For a 0+ld system one has ax = ay such that a± = ax and 'H^t = 'Hzx = 'Hsy, so that one 
obtains 


= 27rNalj dpp^feq(yp^ + 7/3 

= nNal f dp p^ /eq 'H 2 T f ^, 


az m 


X ^xP 

az rh 


'H2,T{OL,m 

UzhioL, rh) = 2TiNal j dpp^feq(yp‘^ + j PC 2 L ^ 


^xP 
az m 


^X ^xP 


.~r 4^2 r + (az m 

= 27 riVa^mM dpp^ - . . ^2 —, 


'0 


+ 




l-Lzm{oL,rh) 

Ussiot, rh) = 271 Nal I dppfeq(yp^ + 7^26 f 


^xP 
m 


^xP 


Also, derivatives of "Ha satisfy 


dTU 

dax 

dUs 

daz 

dUs 

dm 


—fir, 

^X 

-, 

az 

zr^m ■ 

m 


For the isotropic equilibrium case, one has 0 ;^ —)■ 1, A —)■ T, and rh ^ rh, 


eq 


'H3,eq(hleq) = 4:71 Nm] 
^3r,eq {rrieq 

^3m,eq('hleq 


eq 


meq-Ai(?fieq) + 3iF2(meq) 


(rheq) = id3L,eq{rheq) = 47riVmgqiF2(?Aeq) , 


= 47riVmgqiF2(meq) 


(BIT) 

(B18) 

(B19) 

(B20) 

(B21) 


(B22) 

(B23) 

(B24) 


(B25) 

(B26) 

(B27) 
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2. Massless Case 


Taking the massless limit of Eqs. (B6) - (B12) one obtains 


7^3(0:)= \\mT-L^i(y.,m) = QNaxCXy / d(p (B28) 

m-S-O Jq VQ!j_/ 

27r 

ii'ixioL) = lim 'H 3 x{ct,m) = 6Nalay / dcj) cos^ (j)'H 2 T( (B29) 
m-s-O Jq \a±/ 

H3 „(q:) = lim'Hs^fcK, rh) = biVa^^a? [ <70 sin^ 0'H2T f—\ (B30) 

m-s>o y Jq Va_L/ 

= limn3L{ct,m) = 6Naxay c?0ai'H2L f—), (B31) 

m-s>o " Jq Va_i_/ 

'Hzmi.OL) = lim UzmioL, m) = 0 , (B32) 

m—^0 


where 'H 2 , 2 T, 2 L{y) = Td 2 , 2 T, 2 L{y, 0). In the transversally-symmetric case, Ux = ay and T-L^t = 
Td^x = disy, and the functions above simplify to 


nsia) = l2nNalH2 {—), (B33) 

\ax/ 

nsrict) = QTTNal'H2T (—), (B34) 

Kaxd 

'H3l{cx) = UTiNalmi (—). (B35) 

\ax/ 

In the isotropic equilibrium case, one has a* —>■ 1 and A —)■ T, and, as a result, 

?J3,eq = 247riV, (B36) 

’^ST.eq = disLfiq = SvriV . (B37) 
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